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DESIGN OF TWO-DIMENSIONAL CHANNELS WITH PRESCRIBED VELOCITY DISTRIBUTIONS 
ALONG THE CHANNEL WALLIS 1 


By Joun D. SrANrTZ 


SUMMARY 


A general method of design is developed for two-dimensional 
unbranched channels with prescribed velocities as a function of 
arc length along the channel walls. The method is developed 
for both compressible and incompressible, irrotational, non- 
viscous flow and applies to the design of elbows, diffusers, noz- 
zles, and so forth. Two types of compressible flow are con- 
sidered: the general type, with the ratio of specific heats y equal 
to 1.4, for example, and the linearized type, in which y is —1.0. 
Two methods of solution are used: In part I solutions are 
obtained by relaxation methods; in part II solutions are 
obtained by a Green's function. 

Five numerical examples are given in part I including three 
elbow designs with the same prescribed velocity asa function of 
are length along the channel walls but with incompressible, 
linearized compressible, and compressible flow. ТЕ 28' con- 
cluded that if a nonviscous gas with arbitrary y (1.4, for ex- 
ample) were to flow through a channel designed for linearized 
compressible flow (y=—1.0), the resulting velocity distribution 
along the channel walls would be nearly the velocity distribution 
prescribed for the linearized compressible flow. 

One numerical example is presented in part ТІ for an acceler- 
ating elbow with linearized compressible flow, and the time 
required for the solution by a Green’s function in part II was 
considerably less than the time required for the same solution 
by relaxation methods in part I. 


INTRODUCTION 


There are two general types of theoretical problem in 
two-dimensional fluid motion: (1) the direct problem, in 
which the distribution of velocity is determined for a pre- 
scribed shape of boundary, and (2) the inverse problem, in 
which the shape of boundary is determined for a prescribed 
distribution of velocity along the boundary. The direct 
problem is an anlaysis problem; the inverse problem is a 
design problem. This report is concerned with the inverse, 
or design, problem for two-dimensional, irrotational flow in 
unbranched channels with prescribed velocities as a function 
of arc length along the channel walls. 

The design of channels with prescribed velocities is impor- 
tant because: (1) Boundary-layer separation losses can be 
avoided by prescribed velocities that do not decelerate rap- 
idly enough to cause separation, (2) shock losses in com- 


pressible flow and cavitation in incompressible flow can be 
avoided by prescribed velocities that do not exceed certain 
maximum values dictated by these phenomena, and (3) for 
compressible flow the desired flow rate can be assured. by 
prescribed velocities that do not result in “choke flow" 
conditions. 


Several methods of channel design have been developed 
for particular application (refs. 1 and 2, for example). 


reference 1 a design method is developed for accelerating 
elbows in which the velocity increases monotonically along 
the channel walls. The method is developed for incompress- 
ible and linearized (үс--1.0) compressible flow. The veloc- 
ity distribution along the channel walls is not arbitrary and 
the design method applies to elbows only. In reference 2 
a design method is developed for straight, symmetrical chan- 
nels with contracting or expanding walls. The method is 
developed for incompressible flow and the velocities are 
prescribed not as a function of arc length along the channel 
walls but as a function of circle angle in the transformed 
circle plane. A more general design is suggested in reference 3, 
but no attempt is made to develop and apply the method. 
In the present report a general method of design is devel- 
oped for two-dimensional, unbranched channels with pre- 
scribed velocities as a function of arc length along the channel 
walls. The method is developed for both compressible and 
incompressible, irrotational, fonviscous flow and applies to 
the design of elbows, diffusers, nozzles, and so forth. Two 
types of compressible flow are considered: the general type 


, with arbitrary value of y (1.4, for example) and the linearized 


type with y equal to —1.0. Tn general, if the prescribed 
velocity along one channel-wall differs from that along the 
other, the channel turns so that the downstream flow direc- 
tion is different from the upstream direction. This change 
in flow direction cannot be arbitrarily chosen but depends on 
the prescribed velocity distribution along the walls. Equa- 
tions are developed for computing this change in flow direc- 
tion for an arbitrary prescribed velocity distribution with 
incompressible or linearized compressible flow. Two meth- 
ods of solution have been developed for the design method 
and are presented in separate parts of this report. In 
part I solutions are obtained by relaxation methods (ref. 4). 
This method of solution results in complete information 
concerning the distribution of flow conditions throughout the 


1 Bupersedes МАОА TN 2593, “Design of Two-Dimensional Channels with Prescribed Velocity Distributions Along the Channel Walls. I—Relaxation Solutions” by John 
D. Btanitz, 1952, and МАСА TN 2595, “Design of Two-Dimensional Channels with Preseribed Velocity Distributions Along the Ohannel Walls. II—Solntion by Green's Function" by 


John D, Btanitz, 1052. 
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channel and, in addition, can be used to obtain nonlinear 
solutions for compressible flow with arbitrary values of y. 
In part Il solutions are obtained by means of & Green's 
function. This method of solution is limited to incom- 
pressible and linearized (y=— 1.0) compressible flow, but 
the method is more rapid than relaxation methods, provided 
information within the channel is not required. 

The design method reported herein was developed at the 
NACA Lewis laboratory during 1950 and is part of a doctoral 
thesis conducted with the advice of Professor Ascher H. 
Shapiro of the Massachusetts Institute of Technology. 


PART I 


GENERAL THEORY AND SOLUTION BY 
RELAXATION METHODS 


A general method of design is developed for. two- 
dimensional, unbranched channels with prescribed velocities as 
functions of arc length along the channel walls. The method 
is developed for both incompressible and compressible, irro- 
tational, nonviscous flow. Two types of compressible flow 
are considered: the general type with arbitrary value for the 
ratio of specific heats y (1.4, for example), and the linearized 
type іп which y is equal to —1.0. The solutions in part I 
of this report are obtained by relaxation methods and give 
complete information concerning the flow throughout the 
channel. Five numerical examples are given, including three 
elbow designs with the same prescribed velocity as a function 
of arc length along the channel walls but with incompressible, 
linearized compressible, and compressible flow. 


THEORY OF DESIGN METHOD 


The design method is developed for two-dimensional chan-. 


nels with prescribed velocities along the channel walls. The 
prescribed velocity. is arbitrary except that stagnation points 
cannot be prescribed. This exception limits the design 
method to unbranched channels. 


PRELIMINARY CONSIDERATIONS 


Assumptions.—The fluid is assumed to be nonviscous and 
either compressible or incompressible. The flow is assumed 
to be two dimensional and irrotational. 

The assumption of two-dimensional, nonviscous, irrota- 
tional motion limits the design method in practice to channels 
with thin (negligible) boundary layers, such as exist near the 
entrance to the channel or after a rapid acceleration of the 
flow through a contraction in the channel. Even if the 
i boundary layer is thin, the design method is limited to (and 

finds its most useful application for) prescribed "velocity 
distributions that, from boundary-layer theory, do not decel- 
erate fast enough to result in separation of the boundary 
layer, which separation alters the “effective” shape of the 
channel and completely changes the character of the flow. 

In some channels with fully developed turbulent boundary 
layers, the design method might be expected to yield results 
that are satisfactory, although approximate, because for this 
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type of flow the rotational motion occurs primarily in regions 
close tọ the channel walls. In channel walls with thick or 
fully developed laminar boundary layers the design method 
cannot be used, because not only is the rotation of the flow 
important in most of the channel but, if the channel bends, 
important secondary flows develop that are not considered 
by the two-dimensional design method. 

Flow field.—The flow field of the two-dimensional chan- 
nel is considered to lie in the physical zy-plane where z and 
y are Cartesian coordinates expressed as ratios of a charac- 
teristic length equal to the constant channel width down- 
stream at infinity. (All symbols are defined in appendix А.) 

At each point in the channel (fig. 1) the velocity vector has 
a magnitude Q and a direction 0 where Q is the fluid velocity 
expressed as the ratio of a characteristic velocity equal to the 
constant channel velocity downstream at infinity. For con- 
venience, the velocity Q is related to a velocity q by 


4-- да (1) 


where q is the velocity expressed as a ratio of the stagnation 
speed of sound and the subscript d refers to conditions down- 
stream at infinity. 

The flow direction 0 at each point in the channel is measured 
counterclockwise from the positive z-axis. From figure 1 


dx=ds cos 0 (20) 
dy=ds вір 0 (2b) 


where ds is a differential distance in the direction of Q, that 
is, along a streamline. 


% 


FIGURE 1.—Magnitude and ‘direction of velocity at point іп zy-plano, 
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Stream function and velocity potential lIf the condition 
of continuity is satisfied, a stream function V can be defined 


such that 
dy=pQ dn (3) 


where p is the fluid density expressed as the ratio of a charac- 
teristic density equal to the stagnation density and where dn 
is a differential distance measured normal to the direction 
of Q, that is, normal to a streamline. Along a streamline, 


dn is zero во that from equation (3) the stream function y 


is constant. 
IL the condition of irrotational fluid motion is satisfied, 9, 
velocity potential e can be defined such that 


de=Q ds (4) 


Normal to а streamline, ds is zero so that from equation (4) 
the velocity potential o is constant. Thus lines of constant 
e and у are orthogonal in the physical xy-plane. 

Outline of method.—Solutions of two-dimensional flow 
depend on known conditions imposed along the boundaries 
of the problem. In the inverse problem of channel design, 
the geometry of the channel walls in the physical zy-plane is 
unknown. This unknown geometry apparently precludes 
the possibility of solving the problem in the physical plane 
and necessitates the use of some new set of coordinates, that 
is, & transformed plane, in which to solve the problem. 
These new coordinates must be such that the geometric 
boundaries along which the velocities are prescribed are 
known in the transformed plane. It is also desirable, for 
mathematical simplicity, that the coordinate system in the 
transformed plane be orthogonal in the physical plane. A 
set of coordinates that satisfies these requirements is provided 
by e and у, which are orthogonal in the physical zy-plane and 
for which the geometric boundaries are known constant 
values of y in the transformed qj-plane.. The distribution of 
velocity as a function of e along these boundaries of constant 
v is known because, if 


Q= Q(s) 
18 prescribed, equation (4) integrates to give 

e= e(s) , 
from which equations, 

= Q(o) 


The technique of the proposed method of channel design is 
therefore to obtain a differential equation for the distribution 
of velocity in the oj-plane. The velocity distribution ob- 
tained from the solution of this equation is then used to obtain 
the distribution of flow direction, from which distribution the 
channel walls in the physical zy-plane are obtained directly. 
The differential equation for the distribution of velocity in 
the qyj-plane is nonlinear (for compressible flow with y other 
than —1.0) and is solved by numerical methods (relaxation 
methods). 


DIFFERENTIAL EQUATION FOR DISTRIBUTION OF VELOCITY IN 
TRANSFORMED ¢¥-PLANE 


The differential equation for the distribution of velocity in 


the transformed ey-plane is obtained from the equations for - 


continuity and irrotational fluid motion expressed in terms of 
the transformed coordinates დ and v. 

Continuity—The continuity equation expressed in terms 
of დ and y becomes (appendix В): 


1 LR log, Py 0 ов, 2105. ב(‎ (5) 


Irrotational fluid motion.—The equation for irrotational 
fluid motion, expressed in terms of e and y, becomes (ap- 
pendix B): | 

О log, © 00 
OV Әф 


Differential equation for distribution of velocity.—The 
second-order partial differential equation for the distribution 
of log, Q in the transformed ¢)-plane is obtained by differenti- 
ating equations (5) and (6) with respect to e and y, respect- 


076 
ბიბი Thus, 


=0 (6) 


ively, and combining to eliminate —— 


Оз log, p , Ò? ა”. Q o ו‎ e (2-282 log. p | 0 108. d 
F ———.-HL————- + T 


g E] 
р? 9 fog. € dE g ane Pp? LU x 9. 0 (7) 


Equation (7), together with a relation between ი, Q, and qa, 
determines the distribution of log, Q in the ¢y-plane for com- 
pressible flow with a given value of qa and for arbitrarily pre- 
scribed variations in log, Q along the boundaries of constant y. 

Density.—The density p is related to the velocity q by 
(ref. 5, p. 26, for example) 


š _ 1. 
ი=(1-%X– әу” (8а) 


which, from equation (1), becomes 


PES 
ი=(1 Bin Q? אל‎ 
2 


Equation (8b) relates the density ი to the velocity Q for a 
given value of qa. 

Incompressible flow.—For incompressible flow р is con- 
stant and equal to 1.0 so that equation (7) becomes 


(8b) 


2 2 | 
оен о (9) 


Equation (9) determines the distribution of log, Q in the 
оү-рівле for incompressible flow. 


CHANNEL-WALL GEOMETEY 


After equation (7) or (9) has been solved to obtain the 
distribution of log, Q in the transformed ev-plane (for the 
arbitrary specified variations in log, Q with e along the 
boundaries: of constant V), the geometry of the channel 
walls: in the physical zy-plane can be determined from the 
resulting distribution of flow direction 0. 


156 


Flow direction 0.—' The distribution of flow direction 0 
along a streamline (constant y) is obtained from equa- 
tion (6), which integrates to give 


(102) 


where the subscript y indicates that the integration 1s taken 
along a line of constant V and where the constant of integra- 
tion is selected to give a known value of 0 at one value of დ 
along each streamline. The integrand in equation (102) is 
obtained from the distribution of log, Q, which is known 
from the solution of equation (7) or (9). 

The distribution of flow direction 0 along a velocity- 
potential line (constant დ) is obtained from equation (5), 
which integrates to give 


s=- | + (8e log, p 49 XC es) dy 


where the –_ დ indicates that the integration is taken 
along a line of constant o and where the constant of integra- 
tion is selected to give a known value of @ at one value of 
, v along each velocity-potential line. As for equation (10a), 
the integrand in equation (10b) is known from the distribu- 
tion of log, Q obtained from the solution of equation (7) 
or (9). 

Channel-wall coordinates.—The variation in 2 along a 
line of constant V in the ey-plane is given by 


(10b) 


а dz ds 
ds do/, 
which, combined with equations (2a) and (4), integrates to 
give 
cos 0 ` 
= f “0: (118) 
Likewise, 
ain 0 
—— d 11b 
2 r0 V (11b) 
E: 1 Eno (110) 
, Q 7” 
cos 6 : 
ë= d 11 
y s0 4 - (119) 


where the constants of integration are selected to give known 
values of z or y at one value of e along each streamline or at 
one value of № along each velocity-potential line. Equations 
(11а) to (11d) determine the distribution of z and y in the 
transformed ¢gy-plane or, which is the same thing, the shape 
of the streamlines and velocity-potential lines in the physical 
zy-plane. In particular, equations (118) and (110) when in- 
tegrated along the boundaries of constant у in the ey-plane 
determine the shape of the channel walls. 

Turning angle.—In general, if the prescribed velocity dis- 
tribution along one channel wall differs from the distribution 
along the other wall, the channel deflects an amount 49, 
which is the difference in flow direction far downstream and 
far upstream of the region in which the prescribed velocity 
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distribution varies. In part II it is shown. that for incom: 
pressible flow the turning angle A0 is given by 


A0—04— 0, 


LACES), CRD] 


where the subscript 4 refers to conditions upstream at infinity 

-and where the subscripts 0 and 1.0 refer to the channe. 
boundaries along which y equals 0 and 1.0, respectively. A 
similar equation will be given later for the case of linearizec 
compressible flow. 


LINEARIZED COMPRESSIBLE FLOW 


The nonlinear differential equation (7) for the distributior 
of velocity in the gp-plane with compressible flow becomes 
linear and is considerably simplified if a linear variation ir 
pressure with specific volume (1/p) is assumed. ‘This linea: 
relation between pressure and specific volume was first sug- 
gested by Chaplygin (ref. 6) in order to linearize the differ- 
ential equations for two-dimensional compressible flow ir 
the hodograph plane. 

Density.—lf a linear variation in pressure with specific 


= volume is assumed, the density ი” is related to the velocity 


q* by (appendix C) ` 


p* — (1+ q*") —1/2 (13; 
where 
= kp (138. 
and · 
q* =k (13b, 


where the constants k, and k: have been determined so that 
values of p given by equation (13) equal the values of , 
given by equation (88) for any two selected values of q (des- 
ignated by gq and g»). Thus, 


Ее) 
1 psg 
k= ee (140) 


and 


(145 


where p, and p, are determined by equation (8a) for the se- 
lected values of qa and g, respectively. A discussion of the 
selection of qa and qs is given in appendix C. It will be 
noted that, if y is equal to —1.0, equation (8a) has the same 
form. as equation (13). 

Stream function and velocity potential.—For the case o: 
linearized compressible flow it is convenient to define the 
stream function у* and the velocity potential »* by 

dy*=p*g* dn (15. 


ғ 


апа 


de*—=q* ds (16: 
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Continuity.—The continuity equation expressed in terms 
of *ק‎ and y* becomes (appendix D) 


20 А 
азва 2—0 (17) 
where 
q* 
M ee 18 
retro ₪ 
or, conversely, 
ჯ 2u : 
(“та (19) 


Irrotational fluid motion.—The equation for irrotational 
fluid motion, expressed in terms of o* and y*, becomes 
(appendix D) 


Olog,u 00 _ 
S By" = 


Differential equation for distribution of log,u.—The partial 
differential equation for the distribution of log, u in the 
e**-plane is obtained by differentiating equations (17) 
and (20) with respect to e* and y*, სასას and com- 


bining to eliminate 55 7 Thus 


970 
Dp OYE 
o? x ов маа 


Equation (21) determines the distribution of log, v in the 
y*y*-plane for linearized compressible flow with a given 
value of qa and for arbitrarily prescribed variations in 
log, Q, related to log, u by equations (1), (13b), and (18), 
along the boundaries of constant /*. Equation (21) is 
linear and 19, like equation (9) for the case of incompressible 
flow, the equation of Laplace. Thus an incompressible 
flow} solution for the distribution of log, Q in the ey-plane is 
also a linearized compressible flow solution for the distribution 
of log, шіп the y*y*-plane. The transformation from the 
yy-plane is different, however, from the transformation from 
the ¢**-plane so that different channel shapes result in the 
23/-plane. , 

Flow direction 0.—The distribution of flow direction 0 
along а streamline (constant y*) is obtained from equa- 
tion (20), which integrates to give 


0 (21) 


! 


det (228) 


Likewise, the distribution of flow direction 0 along a velocity- 
potential line (constant о%) is obtained from equa- 
tion (17), which integrates to give 


(22b) 


Equations (22а) and (22b) for linearized compressible flow 
correspond to, and are used in the same manner as, equa- 
tions (10a) and (10b) for the usual type of compressible or 
incompressible flow. 


Channel-wall coordinates.—The variation in ₪ along a 
line of constant y* in the e*/*-plane is given by 


oz (ат 7) 
00* de аф“, 


which combined with equations (28) and (16) integrates 
to give 

(23a) 
Likewise, 


(23b) 


(230) 


(23ძ) 


Equations (238) to (28d) determine the distribution of z and 
y in the transformed ¢*y*-plane or, which is the same thing, 
the shape of the streamline and velocity-potential lines in the 
physical zy-plane. In particular, equations (23a) and (280), 


' when integrated along the boundaries of constant V” in the 


e*y*-plane, determine the shape of the channel walls. 
Equations (23a) to (23d) for linearized compressible flow 
correspond to, and are used in the same manner as, equa- 
tions (11а) to (11d) for the usual type of compressible ог 
incompressible flow. 

Turning angle.—In part II it is shown that for linearized 
compressible flow the turning angle, or difference in flow 
direction far downstream and far upstream of the region in 
which the prescribed velocity distribution varies along the 


channel walls, is given by 
2398) |ae* (24) 


мр [Coe ui C oe 


where Ay* is the value of /* along the left boundary (channel 
wall) when faced in the direction of flow if the value of ¥* 
along the right boundary is zero, and where the subscript 
Ay* refers to the boundary along which * is equal to Ay*. 


NUMERICAL PROCEDURE 
'The channel design method in part I of this report was 


developed for three types of fluid flow: (1) compressible, 


(2) incompressible, and (8) linearized compressible. Al- 
though the numerical procedures of the design method are 
similar for each type of fluid, the procedures differ in detail 
and are therefore considered separately in this section. 


COMPRESSIBLE FLOW 


The numerical procedure for channel design with com- 
pressible flow (y=1.4, for example) is as follows: 

(1) The velocity is specified as a function of arc length 
along that portion of the channel walls over which the ve- 
locity varies | 

f 4--4(8) 
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Or qa 18 specified and 


Q=Q(s) | (25) 


where 8 is arbitrarily equal to zero at that point along one 
channel wall where the velocity first begins to vary. 

(2) The channel-wall boundaries of the flow field in the 
transformed ¢gy-plane are straight, parallel lines of constant 
y extending indefinitely far upstream and downstream 
between დ equals + o, where დ is arbitrarily equal to zero 
at that point on the channel wall at which s is equal to zero. 
The value of y along the right channel wall when faced in the 
direction of flow (direction of positive о) is arbitrarily set 
equal to zero in which case the value of y along the left 
channel wall (Ay) is obtained by integrating equation (3) 
across the channel at a position far downstream where flow 
conditions are uniform 

A= pa (26) 

(3) The distribution of log, Q as a function of დ along the 
boundaries m the ev-plane is obtained by integrating equa- 
tion (4) between limits so that 


p= |" Qds=o() 


which together with equation (25) gives the distribution of 
log, Q along the boundaries in the ey-plane 


log, 0=/ (ש)‎ ` (28) 


The integration indicated by equation (27) is carried out 
numerically for arbitrary distributions of Q as a function of s. 

(4) If the velocities prescribed along one channel wall 
differ from those along the other wall, the channel will, in 
general, turn the flow. This turning angle cannot be de- 
termined exactly for compressible flow until the channel 
design is completed. However, it will be shown that this 
turning angle is only slightly greater than that resulting for 
linearized compressible flow with the same prescribed 
velocity and with a suitable selection for qa and g in equa- 
tions (148) and (14b). This latter turning angle for linearized 
compressible flow-is given by equation (24), which can be 
integrated numerically for the arbitrary distribution of 
log, u=f (e) corresponding to equation (28). 

(5) In order to solve equation (7) for the distribution of 
log, Q in the ¢gy-plane, it is convenient to eliminate the 
density terms from equation (7) by means of equation (8b). 
Thus, equation (7) becomes 


Ao log, q v 14 po log, g a 240 (Es log, 9:94) 44D (RA q “о 


(29) 


where 
1 
Ie 
—— —— nu 
7-8 =: 
שר‎ d 


0 


A= 


(27) 
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Equation (29) is nonlinear, and it can be solved by relaxation 
methods (refs. 4 and 7, for example). A grid of equally 
spaced points, at each of which the value of log, Q is to be 
determined, is placed in the flow field between the channel- 
wall boundaries. The grid is extended upstream and down- 
stream sufficiently far so that constant values of log, Q are 
obtained across the channel by the relaxation methods. In 
the numerical examples to be presented six or eight grid 
spaces were used across the channel. In example III the 


number of grid spaces was reduced from eight to four with 


negligible effect on the resulting channel design. The values 
of log, Q at each grid point were relaxed to five significant 
figures. If the same velocity distribution is prescribed along 
both. walls, the channel 18 symmetrical so that the volocity 
distribution in only one half of the channel need be cote 
mined by relaxation methods. 

(6) After log, Q has been determined at each grid point i in 
the ej-plane, the distribution of 0 is determined by equa- 
tions (10a) and (10b), which are integrated numerically. The 
constants of integration in equations (10a) and (10b) are 
determined to give a specified value of 0 at one point in the 
channel (far upstream, for example). The integrands in 
equations (108) and (10b) are determined by numerical 
methods (tables I to VII, ref. 4, for example) from the known 
values of p and log, Q at each of the grid points. If it is 
desired to know the flow direction along the channel-walls 


. only, equation (108) can be solved along the channel-wall 


boundaries /=0 апа ұ--Ау only. If it is desired to know 
0 everywhere in the channel, the recommended procedure 18 to 
determine the variation in 0 along the mean streamline 
(y = (Ay) [2) by equation (10a) and to determine the 
variation in 0 along each velocity-potential line from the 
previously determined values on the mean streamline by 
equation (10b). 

(7) After the distributions of log, Q and 0 are known in the 
ey-plane, the shapes of the streamlines and the velocity- 
potential lines in the physical zy-plane or, which is the same 
thing, the distributions of 2 and y in the transformed yy- 
plane are determined by the numerical integration of equa- 
tions (118) to (11d). The constants of integration in these 
equations are determined so that specified values of z and y 
occur af one point in the flow field. Тһе recommended 
procedure 18 to determine the variation in z and y along the 
mean streamline by equations (11а) and (110) and to deter- 
mine the variation in z ahd y along each velocity-potential 
line for the previously determined values on the mean 
streamline by equations (11b) and (11d). If it is desired to 
know the თ and y coordinates for the channel walls only, 
equations (11a) and (110) can be solved along the channel- 
wall boundaries y=0 and y=Ay only. 
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INCOMPRESSIBLE FLOW 


The numerical procedure for channel design with incom- 


pressible flow (р--1) is similar to that just outlined for 


compressible flow, but with the following differences: 

(1) The velocity is specified as a function of arc length by 
equation (25) alone. 

(2) The value of V along the left channel wall (Ay) 18 equal 
to 1.0 instead of the value given by equation (26). 

(3) The distribution of log, Q as a function of ¢ along the 
channel-wall boundaries in the ey-plane is the same as that 
obtained from equations (25) and (27) and given by equa- 
tion (28). . 

(4) The turning angle A0 of the channel is given by 
equation (12). 

(5) The distribution of log, Q in the gy-plane is obtained 
from the solution of equation (9) by relaxation methods. 


(6) After log, Q has been determined at each grid point . 


between the channel-wall boundaries in the ev-plane, the 
distribution of 0 is determined by equations (10a) and (10b) 
as indicated previously for compressible flow, but with p 
= equal to unity. 

(7) After the distributions of log, Q and @ are known in the 
yy-plane, the shapes of the streamlines and velocity-potential 
lines in the physical zy-plane are determined by equa- 
tions (119) to (11d) as indicated previously for compressible 
flow, but with p equal to unity. 

LINEARIZED COMPRESSIBLE FLOW 


The numerical procedure for channel design with line- 
arized compressible flow (y=—1.0) is similar to that pre- 
viously outlined for compressible flow, but with the following 
differences: | 

(1) The velocity q is specified as a function of arc length 
along the channel walls by q(s) or by 44 and equation (25). 
For each prescribed velocity, there are an infinite number of 
linearized compressible flow solutions depending on the 
selected values of qa and дь in equations (148) and (14b). 
However, for values of qa and ძა within therange of gq prescribed 
along the channel walls (and therefore everywhere in the 
channel), the solutions, that is, channel shapes, probably 
differ only in small detail. The best solution 18 that most 
nearly like the nonlinear compressible solution with arbitrary 
value of y (1.4, for example). In the numerical examples of 
this report it is shown that, if qa and дь are equal to the maxi- 
mum and minimum values of g, a good solution results, at 
least if the ratio of these prescribed velocities is not too large 
(2:1 in the numerical examples). On the other hand, if con- 
tinuity is to be satisfied for a gas with the correct value of y 
(1.4, for example) upstream and downstream of the region 
of the channel in which the prescribed velocities vary, then 
qa Nd gy must equal qu and да. . 

After qa and q, have been selected, the velocity distribution 
q(s) is expressed as g*(s) by equation (13b) where kz 18 given 
by equation (14b) so that 


q* —q* (s) (30) 


The velocity 4% is then expressed as u by equation (18) so 
that 
u=u(s) (31) 


In the particular case where the selected value of q, is equal 
to q», the value of kais given by equation (C4b) in appendix C, 
where the significance of this particular case is also discussed. 

(2) The solution is obtained in the transformed o*y*- 
plane where e* and y* are defined by equations (16) and (15), 
respectively. If the value of /* along the right channel wall 
when faced in the direction of 0” is zero, the value of y* along 
the left wall (Ay*) is obtained by integrating equation (15) 
across the channel at a position far downstream where flow 
conditions are uniform - 


Ay* —ps*qa* ЕКЕУ) 


(3) The distribution of log, v as a function of ф* along the 
channel-wall boundaries in the o*J*-plane is obtained by 
integrating equation (10) between limits similar to those 
discussed previously for compressible flow so that 


e= [at ав (33) 
which together with equation (31) determines the distribu- 
tion of log, u along the channel-wall boundaries in the o*y*- 
plane 


log, u=f (e*) (34) 


(4) The turning angle A0 of the.channel is given by equa- 
tion (24). | 

(5) The distribution of log, u in the »**-plane is obtained 
from the solution of equation (21) by relaxation methods. 

(6) After log, u has been determined at each grid point 
between the channel-wall boundaries in the »*y*-plane, the 
distribution of 0 is determined by equations (22a) and (22b) 
in a manner similar to that outlined previously for compress- 
ible flow.  . | 

(7) After the distributions of log, и and 0 аге known in the 
e*V*-plane, the shapes of the streamlines and the velocity- 
potential lines in the physical zy-plane are determined by 
equations (238) to (23d) in a manner similar to that outlined 
previously for compressible flow. The velocities g* in equa- 
tions (23) are obtained from the known values of v, and the 
densities ი” are given by equation (13). , 


NUMERICAL EXAMPLES 


The channel design method has been applied in part I to 
the five examples listed below: 


Type of channel Type of flow ` 


[m | conversing ection ე סמטקמשת‎  _ 
i 
Тапа compres 


V Elbow Compressible (y=1.4) 


EXAMPLE 1 


4 


The first numerical example is the design of a reducing 
section 1n а straight channel such that the upstream velocity 
is half the downstream velocity. The solution is for incom- 
pressible flow. 


+ QO 


3.0. 


0 L0 2.0 


< 


FIGURE 2,—Prescribed velocity distribution as function of arc length along channel wall for 
examples I, III, IV, and V. Equation (35). 


Prescribed velocity distribution—The prescribed velocity 
as a function of arc length s along both channel walls is 
given by | 
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From equations (35) and (30), log, Q is a known function of o, 
which function 19 plotted in figure 3. 

Results.—' Тһе results of example I are presented in 
figures 4 to 7. | 

In figure 4, lines of constant velocity Q and flow direction 
0 are plotted in the transformed Фұ-ріапе. The flow direc- 
tion 0 is constant and equal to zero along the mean streamline 
)/==0.5( , indicating that the center line of the channel is 
straight. Тһе maximum absolute values of 0 occur along the 
channel walls. The solution is svmmetrical about the mean. 
streamline. "The lines of constant Q and 0 are orthogonal. 


Ф 
1.0 2.0 


Q=0.5 (s <0) 
siat. zh (0 >8 > 3.0) (85) 
V 5:'6 27 So 
Q=1.0 (s = 3.0) 
The prescribed velocity given by equation (35) 18 plotted in 
figure 2. 
Equation (35) together with equation (27) results in 
| e—0.58 | (8<0) 
3 4 
e=stig—jog (OS 8 $3.0) (36) | ₪ 
= -- 38 (8 < 3.0) FIGURE 3.—Prescribed distribution of log, Q as function of e along ond walls for examplo I, 
დ 
-ნ -4 -2 0 .2 4 6 8 lO 12 4а 16 18 20 22 24 2.6 
1.0 AX — 8 ni _ w 7 LO 
< 
.8 | 8 
14 aa [ өз ° 
Q O G 0 2 2) .75) . B 2 2 
y Ep жені на Бей жан = , 
: ЕВЕ , 
= ₪ 
; „იეს. ру ты 
| SEERA - | 
, C SCRHBSEEDRHAIDE A ЛМ А Ж 
-6 -4 -2 0 .2 4 6 8 סו‎ I.2 5 І6 18 20 22 24 2.6 
: | 


FIGURE 4.—Lines of constant velocity Q and flow direction 0 in transformed ey-plane for example 1, Incompressible flow; prescribed velocity given in figure 2. 
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FIGURE 6,—8treamlines and veloolty-potential lines on physical ry-plane for example I. Incompressible flow; prescribed velocity given іп figure 2. 


In figure 5, lines of constant z and y are plotted on the 
transformed ej-plane. Along the mean streamline (ү/--0.5) 
the value of y is constant and equal to zero indicating, as 
before, that the center line of the channel is straight. The 
lines of constant z and y are orthogonal, and the system of 
curves forms a square network. The solution is symmet- 
rical. 

321096—66— 12 : 


In figure 6, lines of constant e and y (velocity potential 
and streamlines, respectively) are plotted in the physical zy- 
plane. The shape of the channel walls is that required to 
result in the prescribed velocity distribution given by equa- 
tion. (35) and plotted in figure 2. The downstream channel 
width is 1.0 by definition. The upstream channel width is 
2.0 in order that the upstream velocity be half the down- 
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FIGURE 7.—Lines of constant velocity Q and flow direction 0 in physical zy-plane for example I. Incompressible flow; prescribed velocity given in figure 2, 


Asymptotic to 
Q= O 


-2 -1 0 l 2 3 4 5 


¢ 


FIGURE 8.—Prescribed volocity distribution as function of arc length along channel wall for 
example П. Equation (38). 


stream velocity. As usual, the streamlines and velocity po- 
tential lines are orthogonal and, with equal increments of 
დ and V, form a square network for incompressible flow. 

In figure 7, lines of constant Q and 6 are plotted in the 
physical zy-plane. The lines of constant Q and 0 are orthog- 
onal. 


EXAMPLE II 


The second numerical example is the design of a converging 
section that funnels the fluid from an infinite expanse into a 


| e= (s— 2) 


straight channel of unit width. Far upstream the channel 
walls are straight and converge at a 90° angle. The solu- 
tion is for incompressible flow. 

Prescribed velocity distribution.— The prescribed velocity 
98 8 function of arc length s along both channel walls is given 
by | 


= (§<0) 
1, 1 177 3 1 /2 
Ое 809 o) a T) 


Q=1.0 


(0 <8 <4)((37) 


(82.4) 


The prescribed velocity given by equation (37) is plotted in 


figure 8. 
Equation (87) together with equation (27) results in 


—2 8 
e=— log. (1-5) (8 <0) 
1 18 1/7 3\83 
7 
| 1 /2 8% (38) 
ag (2-1 / (0 <s <4) 
აეე” (824) 
: дт жы 


From equations (37) and (38), log, Q is а known function of 
,שי‎ which function is plotted in figure 9. 
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«თ 


-0 
FIGURE 9,—Prescribed distribution of log, 0 as function of დ along channel walls for example Н. 


Results.—The results of example 11 are presented in 
figures 10 to 12. š 

In figure 10, lines of constant velocity Q and flow direction 
0 are plotted in the transformed ey-plane. The flow direc- 
tion @ is constant and equal to zero along the mean stream- 
' line )/=0.5( , indicating that the center line of the channel 
is straight. The solution is symmetrical about the mean 
streamline. As for example І, the lines of constant Q and 0 · 
are orthogonal. | 

In figure 11, lines of constant e and v are plotted in the 
physical zy-plane. The shape of the channel walls is that 


-8 -6 -4 se 0 ie 4 6 8 LO 


required to result m the prescribed velocity distribution 
given by equation (37) and plotted in figure 8. As usual, the 
streamlines and velocity-potential lines are orthogonal and, 
for incompressible flow with equal increments of დ and y, 
form a square network. 

In figure 12, lines of constant Q and 0 are plotted in the 
physical zy-plane. The lines of constant Q and 0 are orthog- 
onal. 


EXAMPLE III 


The third numerical example is the design of an elbow for 
which the upstream velocity is half the downstream velocity. 
The prescribed velocities are such that no deceleration occurs 
anywhere along the channel walls. The solution is for in- 
compressible flow. 

Prescribed velocity distribution.—Along both walls up- 
stream of the elbow the velocity Q is equal to 0.5, and along 
both walls downstream of the elbow Q is equal to 1.0. The 
transition from Q equals 0.5 to 1.0 along both walls of the 
elbow will be the prescribed velocity distribution as a func- 
tion of arc length given by equation (35) for example I and 
plotted in figure 2. In terms of log, Q as a function of e, this 
prescribed velocity distribution is given by equation (36) 
and is plotted in figure 3. Although this velocity distribu- 
tion is the same for both walls, the distribution on the outer 
wall (wall with larger radii of curvature) is shifted in the 
positive დ direction an amount equal to 2.25 relative to the 
distribution on the inner wall. 'Thus, & velocity difference 
exists on the two walls at equal values of e, as shown in 
figure 13. The greater this difference in velocity and the 
greater the range in e over which velocity differences exist, 
the greater is the elbow turning angle. For the prescribed 
velocity distribution given in figure 13, the elbow turning 
angle given by equation (12) was 89.37° compared with a 
value of 89.36° obtained from the relaxation solution. 

Results.—The results of example III are presented in 
figures 14 to 16 and im tables I and II. (The numerical 


L2 14 16 18 20 22 24 26 28 2 
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FIGURE 10.—Lines of constant velocity Q and flow direction 0 in transformed ey-plane for example IL Incompressib]e flow; prescribed velocity given In figure R. l 
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FIGURE 11.—8treamlines and velocity-potential lines in physical zy-plano for example II, Incompressible flow; prescribed velocity given in figure 8. 
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. 2 
FIGURE 12,—Lines of constant velocity Q and flow direction 0 in physical zy-plane for example 11, Incompreesible flow; prescribed velocity given in figure 8. 
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Inner woll 
(p and y* 
equal zero) 


FIGURE 13.—Prescribed distribution of log, Q as function of თ along channel walls for examples 
III, IV, and V. š 


results for examples ІШ, IV, and V are tabulated in tables 
I to VI to enable a detailed comparison of the three elbow 
designs with tbe same prescribed velocity Q distribution as 
& function of arc length but with incompressible (example 
ПІ), linearized compressible (example TV), and compressible 
(example V) flow.) 

. In figure 14, lines of constant Q and 0 are plotted in the 
ey-plane. Тһе flow direction 0 varies along the mean 
streamline (/—0.5),-indicating that the channel is curved. 
The solution is unsymmetrical. #5 for examples I and II, 
the lines of constant Q and 0 are orthogonal. 

In figure 15, lines of constant o and у are plotted in the 
physical zy-plane. The shape of the channel walls is that 
required to result in the prescribed velocity distribution 
given by equations (35) and (36) and plotted in figures 2 
and 13. 'The upstream channel width is twice the down- 
stream width in order that the upstream velocity be half 
the downstream velocity. It is interesting to note that, 
before curving in the direction of the elbow turning angle, 
the inner wall first curves in the opposite direction. This 
behavior of the inner-wall geometry is necessary in order to 
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maintain the prescribed constant velocity along the outer 
wall where the velocity would otherwise decelerate because 
of the necessary curvature in the direction of elbow turn- 
ing. This feature of the elbow geometry will also be noted in 
examples IV and V. As usual, the streamlines and velocity- 
potential lines are orthogonal and, for equal increments of 
დ and y, form a square network. 

In figure 16, lines of constant Q and @ are plotted in the 
physical zy-plane. The lines of constant Q and 0 are 
orthogonal. 

EXAMPLE IV 


The fourth numerical example is the design of an elbow 
with the same prescribed velocity Q, as a function of arc 
length, used in example III but for linearized compressible 
flow (y —1.0). 

Prescribed velocity distribution.— The prescribed velocity 
distribution Q is the same 98 that for example III and with 
qa equal to 0.80176. "The variation in Q with s along one 
channel wall is plotted in figure 2. The values of q and თ 
in equations (14а) and (14b) are equal to q, and qa, or 
0.40088 and 0.80176, respectively. For these values of q 
and Ф and for the prescribed velocity distribution with 
linearized compressible flow, the elbow turning angle given 
by equation (24) was 104.08? compared with & value of 
104.07? obtained from the relaxation solution and a value of 
89.36? obtained for incompressible flow (example IIT). 

Results.—The results of example IV are presented in 
figures 17 to 19 and in tables 111 and IV. 

In figure 17, lines of constant 0 and 0 are plotted in the 
transformed e*/*-plane. The solution is unsymmetrical 
and the lines of constant g and 0 are orthogonal. l 

In figure 18, lines of constant y*/Ay* and y*/Ay* are 
plotted in the physical zy-plane (where the constant Ay* 
is given by equation (32) and is equal to 0.73782 for qa 
equal to 0.80176). The shape of the channel walls is that 
required to result in the prescribed velocity distribution used 
in example III but with linearized compressible flow and 
for 04 equal to 0.80176. From continuity considerations 
the upstream channel width is 1.5385 times the downstream 
width. As in example III, the inner wall of the elbow first 


18 20 2.2 24 26 28 30 3.2 34 36 3.8 40 42 44 46 4.8 
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FIGURE 14.—Lines of constant velocity Q and flow direction 6 in transformed хұ-р1апе for example ІП. Incompressible flow; prescribed velocity given in figures 2 and 13. 
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FIGURE 15.—8treamlines and velocity-potential lines in physical zg-plane for example IIl. Incompressible flow; prescribed velocity given in figures 2 and 13, 
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FIGURE 16.—Lines of constant velocity Q and flow direction 8 іп physical zy-plane for example ПІ. Incompressiblo flow; prescribed velocity given in figures 2 and 13. 
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FIGURE 17,—Lines of constant velocity q and flow djrection біп transformed »*y*-plane for example IV. Linearized compressible flow; prescribed velocity as function of aro length along channel 
walls same as for example III (fig. 2) and with 04 equal to 0.80176. 


turns in the opposite direction to the elbow turning angle. 
As usual, the streamlines and velocity-potential lines are 
orthogonal, 


In figure 19, lines of constant q and 0 are plotted in the 


physical zy-plane. The lines of constant g and 0 are not, in 


. general, orthogonal. 


EXAMPLE V 


The fifth numerical example is the design of an elbow with 
the same prescribed velocity Q, as 9, function of arc length, 
used in examples 111 and IV but for compressible flow 
(ү--1.4). : 

Prescribed velocity distribution—The prescribed velocity 
distribution Q is the same as that for examples 111 and IV 
but with იკ equal to 0.79927. The variation in Q with s 
along one channel wall is plotted in figure 2. 

Results.— The results of example V are presented in 
figures 20 and 21 and in tables V and VI. 

In figure 20, lines of constant y/Ay and y/Ay are plotted in 
the physieal zy-plane (where the constant Ay is given by 
equation (26) and is equal to 0.71054 for дг equal to 0.79927). 
The shape of the channel walls is that required to result in 
the prescribed velocity distribution used in examples 111 and 
IV but with compressible flow (ү--1.4) and for qa equal to 
0.79927. The upstream channel width is 1.5412 times the 
downstream width, and the turning angle is 105.31? com- 
pared with 104.07? for linearized compressible flow (example 
IV) and 89.36? for incompressible flow (example ПІ). The 
streamlines and velocity-potential lines are orthogonal. 

The shape of the elbow for compressible flow (example V, 
fig. 20) is nearly the same as the shape of the elbow for 
linearized compressible flow (example IV, fig. 18). There- 
fore, in figure 21 the contours of the walls for both examples 
are compared. Тһе difference in contours is very small and 
it is concluded that, if a nonviscous gas with arbitrary y (1.4, 
for example) were to flow through a channel designed for 
linearized compressible flow (y=— 1.0), the resulting velocity 
distribution along the channel walls would be nearly the 
velocity distribution prescribed for the linearized compress- 
ible flow, at least if the linearized flow were selected (by 
the choice of 44 and თ) so that the densities were equal for 
both types of flow at the maximum and minimum velocities 


and if the ratio of these prescribed velocities is not too large 
(2:1 in the numerical examples). This conclusion is impor- 
tant because the design method for linearized compressible 
flow is considerably faster than the design method for com- 
pressible flow with y other than —1.0. 


PART I 
SOLUTION BY GREEN’S FUNCTION 


In part II a method of solution for the design of two- 
dimensional channels with prescribed velocity distributions 
along the walls is developed by means of the appropriate 
Green’s function. The method applies to incompressible and 
linearized compressible, irrotational flow. One numerical 
example is presented for an accelerating elbow with linear- 
ized compressible flow and with the same prescribed condi- 
tions as example IV of part I. 


METHOD OF SOLUTION 


The method of sdlution by Green’s function is in conjunc- 
tion with a formula derived from Green’s theorem. 


PRELIMINARY CONSIDERATIONS 


Stream function &.—In part II it is convenient to define 
the stream function by Y, where for incompressible flow 


ძV=> dy (39a) 
and for linearized compressible flow (y=—1.0) 
| тау“ 
ы а 


For both types of flow Y varies from zero along the right side 
of the channel, when faced in the direction of flow, to 2/2 
along the left side. l 
Velocity potential Ф.-Іп part II it is convenient to define 
the velocity potential by Ф, where for incompressible flow 


| d$—7 de (402) 
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FIGURE 18.—8treamlines and velocity-potentia! lines in physical zy-p inearizedicom 
; aš for example III (fig. 2) and with{gs equal to 0.80176. 
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FIGURE 19.—Lines of constant velocity g and flow direction 0 in physical ry-plane for example IV. Linearized compressible flow; preseribed velocity as function of aro length along channel 
walis same as for example III (fig. 2) and with gz equal to 0.80176. 
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FIGURE 20.—Btreamlines and те]осіёу-роѓіепНа] lines in physical zy-plane for example V. Com 
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(421.4); prescribed velocity as function of aro length along channel walls samo as 
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FIGURE 21.—Compartson of channel-wall shapes for compressible flow (example V) with y equal to 1.4 and for linearized compressible flow (example IV) for same prescribed velocity as function 


of arc length along channel walls (fig. 2). 
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and for linearized compressible flow 
(40b) 


Channel-wall coordinates.—From part I the distribution 
of channel-wall coordinates დ and y along the boundaries of 
constant Y equal to 0 and я/2 in the transformed $¥-plane is 
given by 


NC ay | раз dé (418) 
Яя ғ 0 
and 
— 2.» oa 
yal ay | ZF ae (41b) 
for linearized compressible flow, and for incompressible flow 
18 given by 
2 Í cos8 
=> I hers ძ% (428) 
and 
_ 2 f sin 0 
וב‎ “0. do (42b) 


where the constants of integration are selected to give known 
(specified) values of z or y at one value of ë along each bound- 
ary. Because g* and Q are known functions of დ from the 
prescribed velocity as & function of arc length along the 
channel walls, the shape of the channel walls in the physical 
zy-plane is given by equation (41) or (42) if 0 is determined 
as a function of Ф along the channel walls. In part II the 
solution for 0 as a function of Ф along the channel walls in the 
@¥-plane is obtained by Green's function. 


SOLUTION BY GREEN’S FUNCTION 


Continuity.—From part I the continuity equation becomes 
in the transformed $Y-plane 


5% tag (43a) 
where for incompressible flow 
V=Q (43b) - 
and for linearized compressible flow 
pP. (430) 
ee | 


Irrotational motion.—From part I the equation for irro- 
tational motion becomes in the transformed S¥-plane 


(44) 


д log. V o0 _ 0 
от ბდ 


Integral equation for 0(d,, Y,).—From equations (438) апа 
(44) 


so that from appendix E the value of 6 at a point (,,¥,) 


ქ 
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within, or on, the channel walls in the transformed &Y-plane 
is given by the integral equation 


009 == I. IC SIT) (0 = | ძთ (46) 


where the subscripts 0 and 2 refer to the channel-wall bound- 


aries along which Y is 0 and 2 respectively, and @ is the 
Green’s function of the second kind for the channel, which 
is an infinite strip of width 7 extending in the $-direction to 


+œ. . 
Green's function G.—The Green's function of the second 
kind 6 for the infinite channel in the ¥-plane is given along 


the channel-wall boundaries (Y equals 0 and z) by (appen- 
dix F) 


Goor == —log, [сов (Ф —Ф,) — cos? (У —Ү,) (47) 


where (Ф, Y) is any point on the channel-wall boundary and 
(P,,) is the point in the channel ог on the boundary]at 
which 0 is to be determined. 


Numerical integration for 6(Ф,7о).-Етот equations (46) 
and (47) 


0 (= | | UA log,V log, [cosh? (6—&,)— 


uns #11 4в-»- 


2 


1. DU : SET log, [cosh? (Ф--%,) - 


cos? #)1 а (ФФ, (48) 


` in which Ње independent variable of integration has been 


changed from dë tó d($—4,) so that the origin, for purposes 
of integration, lies at Ф, rather than Ф--0. If for small 


changes in (Ф--Ф.), that is, for small АФ, the term 2 


may be considered constant and equal to its average a 
over the interval АФ, then 


ò log, V Alog,V 
Ф АФ 


and equation (48) becomes 


თ )5--,(-]-4% 
236 (5, Y) = > [EU log, V | log, [cosh? (6—®,) — 
| (Ф- e - (Ф--Ф.) 


ane რითი». 


2 


T (а log, V ('e—*3+ae 


as NS ტდ 


cos? Vა| d е) 


log, [cosh? (Ф –– დე) — 


($—9,) 


(49) 
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where the summation sign is understood to mean that the 
quantity within the braces is summed over the entire range of 
(დ–– ზე) between + o. 

Equation (49) determines 0 at any point in the flow field 
(channel). For a point (®,,¥,) on the channel walls v, is 
equal to 0 or 7/2 and the integrands in equation (49) become 


2 log, cosh |(@—®,)| 
or 
2 log, sinh |(Ф—%Ф,) |. 


80 that equation (49) becomes 


тоб) > (“> logeV ar wx log, V A (508) 
where 
41 = 1-5 4877 La-a. (50ხ) 
I: aif Y,—0 
1+= 8 if Vა=– 
: 2 (500) 
1,= თ if Y,—; 
where 
I | 
აა + ми нее კინ, cosh (6-2 თIC–თე| (50d) 
ә f 
pot | ор sinh |)9--80|8|66--50| (609 


where the + signs apply for positive values of (6—%,) and 
the — signs apply for negative values of (®—®,). Methods 
of evaluating თ and 8 are given in appendix G, and tabulated 
values are given for a wide range of |)6--60( | in table VII. 
Equation (50a) determines 6(@,,¥,) at any pomt on the 
channel-wall boundaries. Thus from equations (41a) and 
(41b) or (428) and (42b) the coordinates for the channel-wall 
shape in the physical zy-plane can be determined: 


NUMERICAL PROCEDURE 


The numerical procedure for the channel design solution 
by Green's function is the same, except for minor details, 
for incompressible and linearized compressible flow. The 
stepwise procedure is outlined as follows: 

(1) For incompressible flow the velocity Q and for linear- 
ized compressible flow the velocity 0, or which is the same 
thing the velocity Q and the constant downstream velocity 
qa, are specified as functions of arc pen along the channel 


walls 
Q— Q(s) 


4--4(8) 


(51a) 
or 
(51b) 


where s is arbitrarily equal to 0 at that point along one chan- 
nel wall where the velocity first begins to vary. 


(2) Compute V as a function of s from equations (43b) and 
(51a) for incompressible flow or from equations (13b), (14b), 
(430), and (51b) for linearized compressible flow 


V=V (a) (52) 


(3) Compute Ф as a function of 8 from equations (4) and 
(40a) for incompressible flow or from equations (16), (32), 
(40b), and (51b) for linearized compressible flow. In equa- 
tion (32) იკ” is obtained from equations (8a), (138), and (14a). 
For arbitrary distributions of Q or q equation (40a) or (40b) 
is integrated numerically by using, for example, Simpson’s 
one-third rule. Thus 


Ф--Ф(8) | (53) 


(4) From equations (52) and (53) V and Ф are known 
functions of s so that 


V=V(@) (54) 


Thus V is a known function of 6 along the channel-wall 
boundaries in the transformed $Y-plane. 

(5) If the prescribed velocity distribution along one wall 
is different from that along the other, the channel will, in 
general, turn the flow. This turning angle #0 is given by 
equation (H5) in appendix Н. If the turning angle is unsat- 
isfactory, a new distribution of velocity as a function of 8 
(eqs. (51a) and (51b)) is prescribed and steps (1) to (5) 
repeated until the desired value of A9 is obtained. Equation 
(H5) is integrated numerically by using Simpson's one-third 
rule, for example, and equation (54). 

(6) The channel-wall boundaries are straight parallel 
lines of constant V equal to 0 and z/2, and extending to + თ 
in the $-direction. Along these boundaries of constant Y, a 
series of equally spaced points are located at each of which 
the flow direction 0 and the z,y-coordinates of the channel 
walls wil be determined by numerical integration. In 
order to use the tables of თ and 8 presented in this report, 
the point spacing A? must be an even multiple of 7/24. Thus 
the smallest point spacing »/24 is equal to Ха of the channel 
width (z/2). Tor & particular prescribed velocity distribu- 
tion along the channel walls the accuracy of the solution 
increases, and so does the amount of computing, as the 
point spacing 18 reduced. ‘The error for a given point spacing 
depends on the prescribed velocity distribution, and its order 
of magnitude is given by the leading term of the error series 
of the formula used for numerical integration (table VIII, 
ref. 4, for example). For the numerical example presented 
in part II of this report the point spacing АФ was 2/12. 
From equation (54) 


А log, V (log, V)seas— (log, РУ)» (55) 
„ა. АФ š 


where the subscripts Ф and Ф--АФ refer to adjacent points 
along the channel boundaries. 

(7) The value of 0 at each point (%,,%,) on the channel- 
wall boundaries is obtained from equation (50a) in which 
(А log, V)/A® is given by equation (55) and АТ is given by 
equations (50b), (500), and table VII. Note that in equation 


` 


' loge V -l . 
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(50а) the origin has been moved to Ф, by changing from Ф 
to (6—,). Thus the value of (A log, V)/A® for a given value 
of ($—4,) varies with %,. | 

(8) The physical z,y-coordinates at each point on the 
channel-wall boundaries are obtained by the numerical 
integration of equations (42a) and (42b) for incompressible 
flow, or equations (41a) and (41b) for linearized compressible 
flow where Ay* is given by equation (32). The constants of 
integration in equations (41) and (42) are selected to give 
known values of z and y at upstream or downstream positions 
where flow conditions can be considered uniform. 


NUMERICAL EXAMPLE 


The channel design method of part LI has been applied to 
the design of an elbow for the same conditions as example 
IV of part I. The design is for an accelerating elbow with 
no local decelerations of the prescribed velocities along the 
channel walls and with linearized compressible flow. 

‘Prescribed velocity distribution.—The prescribed velocity 
distribution along the channel walls is the same as that for 
example IV of part I. The prescribed velocity as a function 
of Ф is plotted in figure 22. 

Results.—As indicated in table VIII, the elbow design 
resulting from the prescribed velocities given in figure 22 
is the same as that obtained by relaxation methods (fig. 21) 
for the same prescribed conditions (example IV, part I). 

The solution obtained by Green's function (part IL) 
required one experienced computer 3 days, whereas the 
solution by relaxation methods (part I) required about 
10 days. The relaxation solutions provide additional 
information, such as the distribution of velocity across the 
channel; but for the most part this additional information 
is of secondary importance, and the design of channels by 
Green’s function is more rapid and therefore to be preferred 
over the design by relaxation methods. 


-l.2 
-[.3 
-1.4 


FIGURE 22— Varintion іп prescribed values of log, V with Ф along channel walls of numerical 
example іп pert II. 


SUMMARY OF RESULTS AND CONCLUSIONS 


À general method of design is developed for two- 
dimensional unbranched channels with prescribed velocities 
as a function of arc length along the channel walls. The 
method is developed for both compressible and incompress- 
ible, irrotational, nonviscous flow and applies to the design 
of elbows, diffusers, nozzles, and so forth. Two types 
of compressible flow are considered: the general type with 
arbitrary value for the ratio of specific heats y (1.4, for 
example) and the linearized type in which y is equal to 
—1.0. In part I solutions are obtained by relaxation 
methods on a transformed plane the coordinates of which 
are the streamlines and velocity-potential lines in the physi- 
cal plane; in part II solutions are obtained by а Green's 
function. The method of solution in part 1 gives complete 
information concerning the flow throughout the channel, 
whereas the method of solution in part II gives the channel- 
wall coordinates only. 

Five numerical exemples are given in part 1 and the results 
are presented by (1) lines of constant velocity and flow 
direction or lines of constant physical coordinates in the 
transformed plane and (2) streamlines and velocity-potential 
lines or lines of constant velocity and flow direction in the 
physical plane. Among- the fiye -examples are three elbow 
designs for the same prescribed velocity as a function of 
are length along the channel walls but with incompressible, 
linearized compressible, and compressible flow. The numer- 
ical results of these three elbow designs are tabulated to 
enable a detailed comparison of the three designs. 

The shapes of the elbows for compressible flow and for 
linearized compressible flow are very nearly the same; and it 
is concluded that, if a-nonviscous gas with arbitrary y 
(1.4, for example) were to flow through a channel designed 
for linearized compressible flow (y=—1.0), the resulting 
velocity distribution along the channel walls would be 
nearly the velocity distribution prescribed for the linearized 
compressible flow. This conclusion is important because 
the design method for linearized compressible flow is con- 
siderably faster than that for compressible flow. 

One numerical example is presented in part II for an 
accelerating elbow with linearized compressible flow. The 
elbow shape obtained from the solution by Green’s function 
in part IL is the same as that obtained from a solution by 
relaxation methods in part I for the same prescribed con- 
ditions. The time required for the calculations was con- 
siderably less for the solution by Green’s function. 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 


A, B, G, D 
A, B 


TR AN‏ > ס 


coefficients, equation (29) I 

arbitrary constants, equation (Cla) 

Bernoulli’s numbers 

constant, equation (E3) 

Green’s function of the second kind, equa- 
tions (22) and (47) 

integral (თ or B) 

coefficient, equation (14a) 

coefficient, equation (14b) 

length of closed boundary 

distance in zy-plane measured normal to 
direction of flow (expressed as ratio of 
characteristic length equal to channel 
width downstream ₪6 infinity) 

static pressure (expressed as ratio of stag- 
nation density multiplied by stagnation 
speed of sound squared) 

velocity (expressed as ratio of characteristic 
velocity equal to constant channel veloc- 
ity downstream at infinity) 

velocity (expressed as ratio of stagnation 
speed of sound) 

velocity used in linearized compressible 
flow and related to q by equation (13b) 

distance from any point in ®¥-plane to 
point (Ф,, VI.) at which logarithmic singu- 
larity exists 

distance in zy-plane measured along direc- 
tion of flow (expressed as ratio of char- 
acteristic length equal to channel width 
downstream at infinity) 

velocity parameter related to g* by equa- 
tion (18) ы 

velocity parameter defined by equations 
(43b) and (430) for incompressible and 
linearized compressible flow, respectively 

complex functions defined by equations 
(ЕЗ), (Fla), and (F2a), respectively 

Cartesian coordinates in physical plane 
(expressed as ratios of characteristic 
length equal to channel width down- 
stream at infinity) 

complex coordinate, equation (F1b) 

conjugate of z 

integral, equation (50d) 

integral, equation (50e) 

ratio of specific heats 

finite increment 

increment of 


0 


დ and et 


₪0 


Subscripts: 


(Ф--Ф,) 


9 -ך (90--9) 


დ, V, დ", y* , 


әзіз 


flow direction in physical zy-plane (meas- 
ured in counterclockwise direction from 
positive z-axis) 

channel turning angle, equation (12) 

density (expressed as ratio of stagnation 
density) 

density in linearized compressible flow and 
related to ი by equation (139) 

velocity potential used as Cartesian co- 
ordinate in transformed ФУ-різле and 
related to დ or თ” by equation (40a) or 
(40b), respectively 

velocity potential for incompressible and 
linearized compressible flow, respectively, 
equations (4) and (16) 

stream function used as Cartesian coordi- 
nate in transformed ¥-plane and related 
to V or y* by equation (39a) or (39b), 
respectively | 

stream function for Incompressible and 
linearized compressible flow, respectively, 
equations (3) and (15) 

boundary value of ү%, for linearized com- ' 
pressible flow, along left channel wall 
when faced in the direction of flow, 
equation (32) 

any harmonic function in $Y-plane 


quantities related to two velocities (qa and 
Qo, respectively) for which density given 
by equation (88) 1s equal to density p 
given by equations (13), (189), and (13b) 

conditions downstream at infinity 

point in ®¥-plane at which 0 is determined 

conditions upstream at infinity 

left channel wall, when faced in direction of 
flow, along which y* is equal to Ay* 

point at (Ф-Ф) on either channel-wall 
boundary 

point at [(6—®,)+A®] on either channel- 
wall boundary 

along lines of constant о, у, დ”, and у“, re- 
spectively 

right channel wall, when faced in direction 
of flow, along which Y, V, or V” is equal 
to 0 

left channel wall, when faced in direction 
of flow, along which V is equal to 1.0 

Jeft channel wall, when faced in direction 


of flow, along which Y is equal to 2 
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APPENDIX B 


EQUATIONS OF CONTINUITY AND IRROTATIONAL FLUID MOTION IN TERMS OF TRANSFORMED 2, ¥-COORDINATES 


Consider the two-dimensional irrotational motion of a 
fluid particle in the physical zy-plane. The fluid particle is 
defined by adjacent streamlines (constant y) and velocity- 
potential lines (constant о) spaced ёп and ôs apart as 
indicated in figure 23. "The velocity Q is parallel to the 
streamlines and normal to the velocity-potential lines. 

Continuity.—From continuity considerations of the fluid 
particle in figure 23 


0 
> (PQ šn)=0 
עס‎ j 


ბ 6 ბ о 0 L 1 1 oln) 


— Еа s 2 


x 


FIGURE 23.—Fluld particle bounded by streamlines and velocity-potential lines in physical 
rp-plane. 


But, from geometrical considerations (ref. 5, p. 167, for 
example) 


1 06( 08 

п Os Оп (B2a) 
and 

21 (ёв) 980 

9 On 08 (B2b) 


so that equation (B1) becomes 
° logs p 9 Ios e, 2 0р, 106. 4, 5_ 


or 
ბ 5960 р аф, Š o Q de , 00 dy 
7 Әф ТА ში ა dn ძი, _ =e 


which, combined with 800901005 (3) and (4), becomes 
9 log, p 5 +220 ° log, Q 


Equation (5) is the continuity equation expressed in terms 
of o,V-coordinates. 

Irrotaetional fluid motion.—For irrotational motion of the 
fluid particle in-figure 23 · 


= (Q58)— 0 
or 
Q 1 O(88) 
| n às On 0 
But, from equations (B2b) and (B3) 


9 log, Q _ ბმ _ 
БЕЛЕ О8 


მ 210.0 (B3) 


0 


or 
dlog, Q dy 20 ძი _ 


TY —=0 


Op dn до ds 
which, combined with equations (8) and (4), becomes 


0 logs Q _ 90 , (6) 


Equation (0) is the equation for irrotational fluid motion 
expressed in terms of the ¢,¥-coordinates. 


APPENDIX C 
RELATION BETWEEN VELOCITY AND DENSITYZASSUMING LINEAR VARIATION IN PRESSURE WITH SPECIFIC VOLUME 


The approximate, linear relation between pressure р and 
specific volume 1/ი first suggested by Chaplygin (ref. 6) is 
given by 
(Cla) 
from which 
(Cib) 


where A and B are arbitrary constants. 

If р denotes the static pressure expressed as a ratio of the 
stagnation density multiplied by the stagnation speed of 
sound squared, Bernoulli’s equation is 


“P+ qdg=0 
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which combined with equation (Clb) integrates to give the 
approximate relation between velocity and density 


B q 


ӘРІ გ — CODSUMDG (C2) 
For convenience equation (C2) can be written as 
1 *2 
——9g"-—1 
253 4 
ог 
p*—(1 גי‎ )18( 
where 
p*=kip (13a) 
and 
q* = ka (13b) 


The constants k, and k: replace the two arbitrary constants 
in equation (C2), and their values are determined so that 
for any two arbitrary values of თ (designated by qa and თ) 
the values of p given by equation (13) equal the values of p 
given by equation (88). Thus the values of p given by 
equation (13) for 4 equal to qa or ია are correct; for all other 
values of q the values of р are approximate. The constants 
kı and k, are determined from the conditions 


pa — ipa 
qa — oqa 
הי‎ (C8) 
= 
From equation (13) - the conditions given by equation (C3) 


_ NM 


and 


ჩი“ + 
1 s d B b 
“> qv 1- (eet) (146) 
Pods 
where იკ and p, are determined by equation (8a) for the 
selected values of д. and g, respectively. 

The values of qa and 9; might, for example, be selected to 
equal the maximum and minimum values of q (which values 
of q must occur on the channel walls and are therefore known). 
Also, the values of qa and ძა might be selected to equal the 
upstream and downstream velocities 4, and gz. In this case 
the upstream and downstream channel widths would then 
satisfy continuity for a gas with the correct value of y (1.4, 
for example). If the upstream and downstream velocities 
are equal, their value and the value of some other velocity 
(the maximum ог minimum ‘velocity, for example) can be 
selected for qa and ძა; or, if desired, 44 can be equal to თა, in 
which case if 

=q-+e where e—0 


4--4 | 
it can be shown from equations (14а) and (14b) that 


(C4a) 
and 
1 
ks EUIS (C4b) 
2 


This latter case, in which q,—9,—49, corresponds to the 
method used by Chaplygin (ref. 6) and Kármán-Tsien (ref. 8) 


in which the correct relation between p and 2 із replaced 


= ל‎ (14a) by a straight line (eq. (Cla)) that is tangent to the correct 
relation at one point (where q,—9:;). 
APPENDIX D 
EQUATIONS OF CONTINUITY AND IRROTATIONAL FLUID MOTION IN TERMS OF TRANSFORMED დ”, y *-COORDINATES 
Consider the two-dimensional irrotational motion of a О log, ი” dọ“ de" |9 log, g* 22 90 dy" — =, 
fluid particle in the physical zy-plane. The fluid particle is Oe ძვ მთ” dy” ძი 
defined by adjacent streamlines (constant /*( and velocity- 
potential lines (constant თ”) spaced ôn and ôs apart as indi- | OT, from equations (15) and (16) | 
cated in figure 23. The velocity 4% is parallel to the stream- 1 /dlog, p 42 log d 
lines and normal to the velocity-potential lines. 5* ЕЕ ; וו‎ (D1) 
Continuity.—From continuity considerations of the fluid 

particle in figure 23 But, from equation (13) ! გ 

0 91 “ —q dl 

9 (p*g* 5n)—0 1 9log,p* | —q Og, q 
: 08 р“ *ק0‎ Sip” მდ" 
L 

ბ loge p" p* , dlog, g*, 1 oln) | 80 that equation (D1) becomes 
ლგ o ` 98 Ton ბი 08 ვა 
1 % @ y 900 —0 (D2) 

which combined with equation (8929) becomes (144% 9% ' ბ” 


C 
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Finally, if 
T (18) 
u= 
1-414 g* 
then 
| д log, 4 : 
d log. q" ა log, (D3) 
Jg “ს 


so that equation (D2) becomes 

(17) 
Equation (17) is the continuity equation expressed in terms 
of »*, ¥*-coordinates and log, u 


Irrotational fluid motion.—For irrotational motion of 
the fluid particle in figure 23 


Doi er 


Or 


9108, q” 1 C(2§) | --0. 


ES שה‎ Òn 
which combined with equation (B2b) becomes 


d log, q* dy* 90 det 
7 ზს dn მი" 48 


or, from equations (13), (15), and n 


d log, q* 
TP 


Finally, from equations (D3) and (D4) 


О Јов, ш 00 = 
o o 


=0 


779 25 -00 =0 (D4) 
(20) 


Equation (20) is the equation for irrotational fluid motion 
expressed in terms of e*, V*-coordinates and log, и 


-e 


APPENDIX E 


INTEGRAL EQUATION FOR 4@(4,,¥.) 


If the distribution of the angle 0(დ,%) in the transformed 
@¥-plane is harmonic, that is, satisfies equation (45) within 


and on the channel walls (v equals 0 and z) then from 
Green’ s theorem and the theorem of mean value it can be 


shown that the value of 0 at a point (®,,¥,) within (or on) 
the channel walls is given by (ref. 9, p. 204, for example) 


8e 
₪ C ota ы 081) 


where the two integrals on the right side of equation (El) 
represent the line integral around the channel walls in the 
counterclockwise direction with the signs adjusted so that 
2 represents the inner normal to the path of integration. 

The function 6(Ф,7) in ה‎ (E1) is of the form 
(ref. 9, p. 204) 


6(%,%)--2- de 


069 (= 2+ (ა (Фу) (Е2) 


“- 


where r is the distance from any point (%,¥%) to the point 
($,Y,) and where w(%,¥) is an arbitrary function that is 
harmonic within and on the channel walls. (Thus from 
equation (E2), G(®,¥) is harmonic within and on the channel 
walls except at the point (@,,%,) where a logarithmic singu- 
larity exists.) Because the harmonic function ₪)6,( is 
arbitrary, the function G(#,¥) can be selected so that along 


the channel-wall boundaries 6 equals 0 and т) 06 is ₪ con- 


stant 6 given by the following equation (obtained from notes 
presented by Tamarkin and Feller in the 1941 Summer 
Session for Advanced Instruction and Research in Mechanics 
at Brown Univ.): 


с=27 (Ез) 


where 7 is the length of the path along which the line integral 
is taken. For the path under consideration Г is infinite and 
therefore G($,Y) can be selected so that Ее 
channel walls. А function with this property is a Green's 
function of the second kind. Equation (#1) becomes 


I. fe 900 00 
En t= Ін IC S) (? 26 | dà 
, 2 . 
or, combined with equation (43a) 


0 (,, = I. 6 — -)5 (| dé (40) 


Along the channel walls 20:14 is known from the prescribed 
velocity distribution so that, after the proper Green's func- 
tion. G has been determined (appendix F), equation (46) de- 
termines the value of 0 at any point (%,,¥,). The value of 
0($,,Y,) given by equation (46) can be adjusted by an arbi- 
trary constant of integration to give a specified value of 0 at 
one point in the flow field. 


is zero along the 
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APPENDIX F 


GREEN'S FUNCTION OF SECOND KIND 


From appendix E Green's function of the second kind Q 
satisfies the condition 


along the channel walls, which are straight and parallel 
boundaries (v equals 0 and 5) extending to + o'in the Ф- 


direction, and satisfies the equation 


06 о 
зы tow? 


everywhere in the channel except at the point (Ф,,Ү,) 
where G has a logarithmic pole. For these conditions the 
Green's function G can be obtained by analogy from the 
velocity potential for incompressible flow into a point sink at 
(®,,V,) between straight parallel boundaries at T equal to 0 


and = ‘The logarithmic pole for G at (®,,¥,) corresponds 


2 
to the point sink, and the condition 280 at the boundaries 
corresponds to zero velocity, that is, no flow normal to the 
boundaries. 

The velocity potential for fluid flow with the boundary 
conditions just described is obtained from two infinite 
series of point sinks with the sinks of each series spaced » 
distance apart in the Y-direction,and the two series arranged 
by the method of images in such а manner that no flow crosses 
the boundaries, that is, 96 о, This arrangement of point 
sinks is shown in figure 24. 

The complex function w, for the first infinite series of 
point sinks is given by (ref. 10, p. 112, for example) 


w= —log, sinh (z— Zo) (Ela) 


(F1b) 


The complex function 4 for the second infinite series of 


where 
z=+1Ť 


nds (mirror image of thea Aust Acris in order to provee 


flow across the boundaries Y equals 0 and 2) is given by 


w= —log, sinh )2--20( (E22) 


(F2b) 


The complex function 4» for the combined flow becomes from 
equations (Fla) to (F2b) | 


where 
2--Ф-4” 


w=w,-+- —log, sinh (Ф-ӘД--4(%-%ДІ-- 


log, sinh [($— დ,) +i(¥+W,)] (F8) 


Point sinks given by 
iv o wj *-loge sinh |(თ-რე) + у) 
а Wo *-log, sinh [ Ф-Ф,) + ҰЛ) 


т/2 


- Z 
Channel wolls-< x – 


` 
` 
` O 


-т/2 


ст = 
FIGURE 24.—T wo infinite series of point sinks required in the development of Green’s function 


of the second kind 6. 


Тһе Green's function of the second kind G corresponds to 
the velocity potential for the incompressible flow and is there- 
fore given by the real part of equation (ЕЗ) 


G— -- log, [cosh? ($—4;) —cos* (Y—Y.)][cosh? (Ф-Ф)- 


cos? (Yr +T) (F4) 
But along the channel walls Y is equal to 0 or; so that 
сов? (Y ש-1‎ ( = cos? (Y — დე 
and equation (F4) becomes 
боо — 100, [cosh? (Ф —Ф,) — сов? (v —Y,) (47) 


Equation (47) gives the Green’s function of the second kind 
along the channel walls (straight parallel lines of constant Y 


equal to 0 and 2 and extending to + o in the @-direction). 
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APPENDIX G 


EVALUATION OF < AND 8 


Several techniques, depending on the magnitude of the 
upper limit |($—9,)|, were used to evaluate the integrals 
თ and В given by equations (50d) and (50e). 
is treated separately in this appendix, and the values of 
(Ф--Ф,) for the upper limit |($— $,)| are considered positive. 
For negative values of (Ф--Ф,) the magnitudes of J (that 
is, of თ or В) are equal for corresponding values of |(6—,)| 
but opposite in sign. As a result the values of AT have the 


same sign. 
INTEGRAL a 


Small and medium values of (Ф--Ф)).-Еог small and 
medium values of the upper limit of integration (Ф--Ф,) 
in equation (50d), that is, for 0 ($—4&,) <60я/24, the 
integral « is evaluated by Simpson's one-third rule using 
increments of ($—4,) equal 00 6 

Large values of (Ф--Ф,).-Еог large values of (Ф—%,), 
that is, for )6--80(<607/24, the integrand in equa- 
tion (50d) becomes 


log, cosh (®—®,) = )6--90( —log, 2 


so that equation (50d) becomes 
60x / 24 ` 
₪ f log, cosh ($—%; d(#—ë,)+- 


(Ф-% 
წირი 2—2), 2] dp 


64+ 


–_–_ )—25.398552 | 


= 0.411980--0.6981476-- (+5 ©, (G2) 
Equation (G2) gives values of a for values of (Ф--Ф,) 
equal to or greater than 607/24. Values of the integral а are 
tabulated in table VII for a range of |($—4,)| between 0 
and 1007/24 in increments of 4/24. For negative values 
of ($—4,) the sign of o is negative. 


INTEGRAL 8 


Small values of ($-—4d,).—For (Ф-Ф,) equal to zero the 
integrand of equation (50e) becomes infinite so that Simpson's 
one-third rule cannot be used to evaluate B in this region of 
(Ф-Ф,), as was done for a. However, equation (50e) 
integrates by parts to give 


Each integral 


(G1) 


(ФФ 

/ ” log, sinh (C-–-დე 8)8--,(=-)9-8/( log, sinh (6—6,)— 
(Ф2--%,) 

ji (©—&,) ctnh (ФФ) 4Ф-%) (G3) 

where the integrand (Ф— Ф,) ctnh (6—®,) on the right side 

of equation (G3) can be expanded in the following series form: 


€ 3 — d M 
(©—%,) cmh @—s)= 14-2456 A _ P 5 ში L 


258, (თ-- ו‎ დე" 218. (C-- დე. 
———— ი 


213 Bi (b –– დ,)!? 
— ^ ma dl 


where 8), Bs, and. so forth, are Bernoulli's numbers (ref. 11, 
p. 90, for example). From equations (G3) and (G4) 


ЕЕ 2) log, sinh ($—4,)— ($ —4$,)— @— - $—4) _ 


(G4) 


(Ф--%) .( 
“ო ד‎ 555: 


2(6—6,)" G—29' 2(5—85)" — - დ)!პ 
6615 | 42,525 1,029,105 ' §,300,667,375 | 


Equation (G5) was used to obtain В as a function of ($— თ) 
for 0<(—6,) <8x/24. 

Medium values of (Ф-Ф,).--Ког medium values of the 
upper limit of integration (8—%@,) in equation (50е), that is, 
for 8т/24<<(Ф--Ф,) X 602/24, the integral В is evaluated by 
Simpson's one-third rule as was done for a. 

Large values of (Ф — ,).—For large values of (Ф--Ф,), that 
is, for (თ-- რე >60»/24, the integrand in equation (50e) 
Becomes 


. (G5) 


log, sith (+ë) = (6—6,) —log, 2 
so that equation (50e) becomes 


is E Toes inh ey eo Sd 
/ -0 (@—&,)—log,.2] d @—&,) 
60a /24 


E XY. 
~24.576082-+| =F —0.693 147 (Ф—Ф,) — 25.3 98562 | 


(G6) 


~ —0.822470—0.693 147 65-55 (დ-- დე? (G7) 
Equation (G7) gives values of В for values of (b—®,) equal to 
or greater than 607/24. Values of the integral В are tabulated 
in table VII for a range of |(Ф-Ф,)| between 0 and 1007/24 
in increments of 7/24. For negative values of (Ф— Ф,), the 


sign of 8 changes. 


APPENDIX H 
CHANNEL TURNING ANGLE | 


If the prescribed velocity distribution along one channel 
wall differs from the distribution along the other wall, then 
in general the channel deflects an amount A6, which is the 
difference in flow direction far downstream and far upstream 


of the region in which the prescribed velocity distribution 
varies. Thus, 


A0—04— 06, (H1) 
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For large values of |(—®,)| such as occur far upstream | ` 1 =)" ° log, V log, კ 
and far downstream of the region in which the prescribed E o| (2408 .-( ОФ 5( | P (H3) 


velocity varies along the channel walls 
cosh? (დ-–-თნე) > > cos? (¥—Y,) 


80 that from equation (47) 
/ 


Likewise, far downstream თე >თ so that 
|)8--%0( | = (b —,) 


/ = Ge — —2[|(b—4,)|—log, 2] (H2) and equation წი substituted into equation (46) gives 
Far upstream Ф,<Ф so that | 0 p ° log, x 2m 0 log, oe | dà (H) 
|(f—%,)|=(—®,) | 
and because V is harmonic From equations (H1), (H3), and (H4) 


ІШЕ: d log. V | הרה.‎ T mE 2 (= IC ა 9106, V) ] ж 6 


so that equation (H2) substituted into equation (46) gives Equation (H5) determines the channel turning angle A6. 
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TABLE I—DISTRIBUTION OF VELOCITY 0 AND FLOW DIRECTION 6 IN TRANSFORMED ¢}-PLANE FOR EXAMPLE ІП 


(ELBOW WITH INCOMPRESSIBLE FLOW) 
[Prescribed variation in Q with aro length s along channel walls plotted in fig. 2; Qu=0.5, 0კ=1.0, A9=89.36°] 
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TABLE II—DISTRIBUTION OF PHYSICAL COORDINATES z AND y IN TRANSFORMED gy-PLANE FOR EXAMPLE III (ELBOW 


WITH INCOMPRESSIBLE FLOW) 
{Prescribed variation in Q with aro length s along channel walls plotted іп fig. 2; Qu»«0.5, 04-21.0, A09 89.30?] 
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TABLE III—DISTRIBUTION OF VELOCITY q AND FLOW DIRECTION 0 IN TRANSFORMED *yp*-PLANE 
FOR EXAMPLE IV (ELBOW WITH LINEARIZED COMPRESSIBLE FLOW) 


[Prescribed variation in Q with arc length s along channel walls plotted in fig. 2; 0.=0.5, Qe=1.0, qa=0.80176, ტ#” 0.73782, А0=104.07°] 
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TABLE IV—DISTRIBUTION OF PHYSICAL COORDINATES с AND y IN TRANSFORMED орж РАМЕ FOR EXAMPLE IV 


1.0 


М 


(ELBOW WITH LINERARIZED COMPRESSIBLE FLOW) 
[Prescribed variation іп Q with are length s along channel walls plotted in fig. 2; Q.=0.5, 0კ=1.0, ga =0.80176, Ay*=0.73782, A= 104.07°] 
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TABLE V—DISTRIBUTION OF VELOCITY q AND FLOW DIRECTION 0 IN TRANSFORMED «eV PLANE FOR EXAMPLE V (ELBOW 
WITH COMPRESSIBLE FLOW (7=1.4)) . 


[Prescribed variation In Q with arc length s along channel walls plotted in fig. 2; Qa=0.5, Cკ=1.0, q4 0.79927, Ay 0.71054, Ag=105.31°] 
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42/6 . 7993 —105. 31 . 7992 —105. 31 . 7002 —105. 31 . T901 —105. 31 . 7992 —105. 31 „7992 —105. 31 . 7993 —106. 31 
43/0 . 1093 —103. 31 . 7083 —]05. 31 „7992 —105.31 | + .7992 —105. 1 . 7992 —105. 31 . 7993 105 31 . 7903 1 
44/6 . 7993 —105. 31 . 7993 —105. 1 . 7993 —105. 31 . 7993 —105. 31 . 7993 —]05. 31 . 7093 —105. 31 . 7093 —105. 31 
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TABLE VII—TABULATED VALUES OF THE INTEGRALS > AND в FOR A RANGE OF |(#—®,)| 
{Computational methods given m appendix а] - 


N АТ-Ас АГ AB 
102—Ф.)] 
А8 А8 
(Abm ж [24) ) 2x[24) 
0 
1(X/24) —(. 611680 
- 2(#/24) 
3(x/24) —. 242701 
4(x/24) — 
5(ж/24) —. 004079 
~ 6(x/24) š - 
7)]24( . 012092 
8 (x/24) 
9)+]24( . 100549 
10(x/24) 
11(x/21) © . 180181 
12(«/24) 
13(х/24) .255015 
14(x/24) 
i 15(X/84) . 327305 
.398006 . 
„467817 
г 537106 
„ 606089 
. 674890 
„ 743685 
„812215 
. 880808 
940380 
1. 017938 
1. 088488 
. 39(x/24) 1. 155034 А 
40(«/04) 
41(w/24) 1. 223576 
42(v/24) 
> 43(v]21) 1. 292118 
44(x/24) 
45(x/24) 14677370 |---------| - 1. 360858 
48)%]24( 14. 366068 ` - 
47 (x[24) 15. 072101 : 1 429106 
48(ұ/24) 15. 795268 
40(x/24) 16. 535571 1. 497738 


. 757436 


(0 For negative values of (Ф--2,) the signs of > and f change, but the signs of Aa and Af remain unchanged, 
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TABLE VII—TABULATED VALUES OF THE INTEGRALS < AND 8 FOR A RANGE 05 )])-.) 
[Computation methods given in appendix 0.| 


АТ = А8 
| (—*,)| 
(Ate 4) (АФ АЗ np 
50(«/23) 
51(х/24) 
52 (x/24) 
53 (x/24) 
64(#/24) 
55(=/24) 
56(x/24) 
57)+/24( 
53(x[21) 
59(я/24) 
60(»/24) 
61 (ж/24) 
62(x/24) 
63(x/24) 
64 (7/24) 
65)*]24( 
66(=/24) 
67(x/24) 
68 (w/24) 
69 (x/24) 
70(x/21) 
712/24) 
72(х/24) 
T3(x[24) 
74(>/24) 
75)*]24( 
76(ж/24) 
TI (x24) 
78(«/24) 
79(«/24) 
80(х/24) 
81)+/24( 
82)%]24( 
83(x/24) 
84(x/24) 
85)+]24( 
86)+]24( 
87)%]24( 
88(x/24) 
89(»/24) 
90 (17724) 
91(x/24) 
92)+/24( 
93(X/24) 
94 (x/24) 
95(x/21) 
96(x/24) 
97 (x/24) 
88(x/24) 
99 (x/24) 
100(x/24) (2 


THE 


ЛЕ 


ЕН 


| ELE 
| ₪ | ₪ 
3 
g 


p 
. 
მ 
= 
a 


ШІ 


i 
: 


ІШІП 


77. 011599 76. 777899 


€ (D For Нче a огрее cn sia ola витр change, but tho Meis of Ae ond გი roman unghanged. 
@ For values of |(დ--დ,)| >100C-/21) use equation ((12) for a and equation (G7) for p. 
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TABLE VIII—COMPARISON OF ELBOW DESIGNS OBTAINED FROM SOLUTIONS BY RELAXATION METHODS AND BY 
GREEN’S FUNCTION | 


[Linearized compressible flow; prescribed velocity distribution given in figs. 2 and 22.) 


¥=( (Inner wall) M=/ (Outer wall) j 
Ф Solution by relaxation Solution by Green’s function Solution by relaxation Bolution by Green's function 
Q methods (Part I) (Part ID methods (Part 1) (Part II 
q 
0 ნ გ მ ტ ნ 
ә) , (deg) Е d (deg) 2 , (deg) გ (deg) 
0. 4009 —2. 466 —0. 760 0 —2. 466 —0. 769 0 0. 5000 0. 4009 —2. 465 0. 770 0 —2. 466 0 
. 4009 —2 241 —. 769 .01 | —2. 241 —. 769 . 01 . 5000 . 4009 —2.241 . 770 -.0 | —2, 241 —.01 
. . 4009 —2.016 —. 760 .0 | 016 —. 769 MI)! 5000 . 4009 —2. 016 . TIO --.01 | —2.016 --,01 
.4009 —1. 791 -.7 .02 | —1.791 --. 769 . 02 . 5000 . 4009 —]. 791 . 7/8 —. 02 | —1.791 —, 02 
. 4009 —1. 566 —. 768 .03 | —1. 566 —. 768 .03 . 5000 . 4009 - 556 . 770 -.0 | –1.§66 --. 03 
4009 —1. 341 –,.? .05 | —1,341 —. 168 . 05 . 5000 . 4009 —1. 341 . 760 —.05 | —1. 341 —. 05 
4009 —1.116 —. 768 „08 | —1.116 —. 768 . 08 . 5000 . 4009 —1.116 . 769 -.08 | —1.116 —. 08 
„4009 —. 891 —. 768 . 14 —. 891 —. 768 .14 . 5000 . 4009 —. 801 . 769 --.13 —. 891 -.8 
. 4009 —. 666 —. 767 .21 —. 666 —. 767 . 24 . 5000 . 4009 —. 666 „768 —.21 —. 660 —, 22 
. 4009 —.441 -.1 .40 --.441 —. 766 41 . 5000 . 4009 --.441 . 767 —. 35 —. 441 —, 30 
. 4009 —.216 —. 763 . 70 —.216 —. 763 . 74 . 5000 ` .4009 —. 216 . 765 —. 56 —, 216 —, 59 
. 4009 - 008 —. 760 1.31 . 010 —. 759 1. 52 . 5000 . 4009 . 009 . 763 --. 92 . 009 ' -1 
. 4072 . 233 —. 782 2.82 - 233 —. 761 2.76 . 5000 . 4009 . 224 . 758 - 5 .4 —1.48 
. 4243 . 450 —. 739 3. 88 . 449 —. 738 3.76 - §000 „4009 . 459 . 751 --2.22 . 469 <2. 21 
. 4489 . 656 —. d 4.00 „658 —. 724 3. 89 . 5000 . 4009 . 684 . 740 — 3. 28 . 684 — 8, 29 
„4780 . 851 —. 712 3.17 . 850 --.713 2.07 . 5000 „4009 . 908 +725 —4. 65 . 908 —4, 08 
. 5004 1. 033 —. 704 1.44 1. 033 —. 705 1.38 - 5000 . 4009 1.132 . 703 —6. 41 1.132 —(. 43 
. 5415 1. 205 —. 703 —. 98 1. 205 —. 704 --1.04 . 5000 .4009 1.355 . 074 –8. 52 1.355 — 8, 56 
. 5782 1.366 —. 710 —4. 00 1. 307 —. 711 —4. 01 . 5000 4009 1. 577 . 636 - 4 1.577 - 11.00 
. 6038 1. 519 —. 725 --7.48 1.519 --.727 - 51 - 5000 4009 1. 707 . 587 —18. 91 1.797 - 41 
„6329 1. 663 —. 749 —11.35 1. 663 —. 750 —11.37 . 5000 4009 2. 013 . 527 --17.18 2.013 --17.16 
„6602 1. 298 —. 781 —15. 52 1. 799 —. 783 —15. 54 . 5000 4009 2. 226 „455 --20. 67 2. 226 — 20, 69 
. 6855 1. 928 —. 822 —19. 98 1.926 —. 823 —19. 98 . 5000 4009 2. 434 . 368 —24. 40 2.434 — 24, 52 
. 7086 2. 016 —. 871 — 24. 62 2. 016 —. 871 —H. 63 . 5000 . 4009 2. 035 . 208 —28. 2. 635 — 28. 59 
. 7293 2.157 —. 928 —29. 46 2. 158 —. 929 —29.47 . 5000 . 4009 2. 829 1853 —32. 89 2. 828 —32. 00 ` 
.7477 2. 261 —. 993 — 34.45 2. 261 –.99+ --34.46 „ 5000 . 4009 3.013 . 024 —37. 40 3. 012 -1 
. 7636 2358 | —1.001 --29. 56 2356 | —1. 066 —39. 57 . 5000 . 4009 3. 186 —.120 —42. 08 3. 185 —42. 09 
. 7109 2. 443 - 143 —44. 76 2. 443 --1.144 --44.77 . 5000 .4009 3. 346 —. 278 —46. 93 3.346 - 41 
. 7875 2. 521 —1.228 — 50. 02 2 521 —1. 229 - 1 - 5000 - 4009 3.492 —. 449 —51. 93 3.492 - 41 
. 7953 2. 590 —1.318 — 55. 27 2. 590 —1. 320 — 55. 27 . 5000 . 4009 3. 623 —. 032 —57. 08 3. 623 — 0 
. 8001 2. 650 —1.414 –60.49 2. 650 —]. 415 — 600. 47 - 5000 . 4009 3. 730 —. 826 — 62. 44 3. 730 - 9 
. 5018 2. 701 —1. 514 -65 55 2.701 - 516 —65. 1 . 5000 . 4009 3.8840 | —1.030 — 608. 16 3. 829 - 7 
. 8018 2. 743 —1. 619 —70. 32 2 744 —]. 620 - 29 . 5070 .4072 3. 901 —]. 242 —74. 80 3. 901 —74. 15 
. 8018 2. 777 —1. 726 - 41 2777 | —1 727 —74. 78 . 5293 . 4243 3.917 | —1.455 —81. 03 3. 945 --80, 01 
. 8018 2. 803 —1. 836 —18. 97 2. 803 —1.837 - 6 . 5599 . 4489 3. 969 —1. 661 —86. 33 3. 969 -% 
. 8018 2.820 | —1. H7 — 82. 82 2. 821 -—]1.948 —82. 80 . 5962 . 4780 3. 074 —1. 856 — 90, 69 3. 975 — 90, 50 
- 8018 2. $31 —2. 059 — 88. 28 2. 531 —2. 059 — 88. 26 . 6354 . 5091 3.966 | —2.038 —$4. 16 3. 007 —91. 09 
. 8018 2 835 —2.171 —89. 37 2. 836 —2.172 --89. 4 . 6754 . 5415 8.050 | —2.209 — 96. 93 3. 051 -8 
. 8018 2.834 | 283 --92. 07 2. 834 <2. 285 —92. 01 . 7149 . 5732 3.927 | —2.369 —99.11 З. 928 —09. 
. 8018 2. 827 —2. 396 — 1H. 40 2. 828 —2.397 —94. 37 .1 . 6038 3.900 | - 520 —100. 83 3. 901 —100. §0 
. 8018 2. 817 —2 508 —96. 39 2. 817 —2. 509 --96. 36 . 7894 . 6329 3. 871 —2. 603 —102. 17 3. 872 —102. 14 
. 8018 2. 803 --2. 619 —898. 05 2. 803 —2. 620 - 3 „8235 . 6602 8. 841 —2, 799 —103. 10 3, 812 —103. 17 
. 8018 2780 | —2.731 —99. 44 2. 786 —2. 732 —99. 42 . 8550 „ 6855 3. 809 —2. 929 --103. 95 8, 810 —103. 91 
. 8018 2.768 | 811 —100. 56 2.767 | —2.812 | —100.55 . 8838 7085 3.777 | —3.055 | —101.49 3.778 --101,48 
.8018 274 | -2.052 | —101.47 2.745 —2. 953 —101.45 9097 7293 3.746 | —3.176 —104. 81 8, 746 ~101. 83 
„5018 2. 72L —3. 082 | —102.18 2.722 | —3. 083 —102. 17 . 9326 7477 3. 714 - 291 --105. 03 3.715 —105. 03 
. 5018 2. 697 —3.172 | —102.73 2698 | —3.173 —102. 72 0521 7636 3. 633 —3. 409 —105. 10 3. 084 —105. 09 
. 8018 2.672 | —3.281 —103. 14 2. 672 —3. 283 —103. 14 9690 7769 3. 653 - 522 —105. 05 3. 654 —105, 01 
. 5018 2.016 | —3. 391 —103. 45 2. 617 —38.392 | —103.44 0822 7875 3.623 | –3.633 —104. 91 3. 021 —101.91 
. 8018 2. 620 —3.500 | –103.66 2. 620 - 501 —105. 66 9919 3. 501 —3. 744 —101.71 3. 505 —104. 72 
. 8018 2.5093 | —3.609 | —103, 81 2. 591 —3. 611 —103. 79 9979 3. 585 | —3.853 | —104 49 3. 566 —101.40 
.8018 2. 566 —3. 719 --103. 90 2 567 —3. 720 -103.90 | 0 8018 3. 537 —3. 962 —101. 28 3. 538 - 0 
. 8018 2. 539 - 828 —103. 97 2. 510 —3. 829 —103.97 | 0 8018 3. 510 „071 —101.18 3. 510 —101. 19 


